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We analyze the interaction between Dark Energy and Dark Matter from a thermodynamical 
perspective. By assuming they have different temperatures, we study the possibility of occurring a 
decay from Dark Matter into Dark Energy, characterized by a negative parameter Q. We find that, 
if at least one of the fluids has non vanishing chemical potential, for instance fj.^ < and fidm = 
or fix = and fidm > 0, the decay is possible, where fi^ and fidm are the chemical potentials of 
Dark Energy and Dark Matter, respectively. Using recent cosmological data, we find that, for a 
fairly simple interaction, the Dark Matter decay is favored with a probability of ~ 93% over the 
Dark Energy decay. This result comes from a likelihood analysis where only background evolution 
has been considered. 

PACS numbers: 95.36.+x,98.80.Es 



I. INTRODUCTION 

There is no doubt that two of the main challenges of 
the present cosmology concerns the understanding of the 
nature of the Dark Energy (DE) and Dark Matter (DM). 
The former is responsible for the present accelerated ex- 
pansion of the universe, as indicated by type la Super- 
novae observations and represents about 70% of 
its material content. The latter one acts exactly like the 
ordinary matter, but does not interact with it, except 
gravitationally, corresponding to about 25% of the ma- 
terial content of the Universe (see Q for a review). In 
a series of recent papers the possibility of a coupling be- 
tween DM and DE has been considered, and we have 
three possibilities: DM decaying to DE BMSlDE de- 



caying to DM 
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interaction in both directions [18|, ll9| . Interaction among 
these completely different fluids has some important con- 
sequences, like addressing the coincidence problem, for 
instance. The coincidence problem can be solved or al- 
leviated on these models by assuming that the DE de- 
cays into DM, thus diminishing the difference between 
the densities of the two components through the evolu- 
tion of the Universe. Thus, if one finds that it is the DM 
that decays into DE, obviously these models can not be 
used to alleviate the coincidence problem, but, on the 
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other hand, it is the DM that dominates in the past, and 
as the decay occurs, the dark energy starts to dominate, 
as indicated by observational data. Thus, the possibility 
of interaction between DE and DM can be investigated 
and the final answer shall be given by the cosmological 
observations. 

Although the observational constraints indicate the 
possibility that the decay can be in both directions, the 
last mechanism, namely the decay of DE into DM, has 
been only recently studied from a thermodynam- 

ical point of view and it seems to be strongly favored by 
the second law, provided that the chemical potential of 
both components vanish. In this paper wc show that, if 
the chemical potential of at least one of the fluids is not 
zero, the decay can occur from the DM to DE, with no 
violation of the second law. The inclusion of a non null 
chemical potential to the dark energy fluid is not a merely 
theoretical tool. As showed in some recent papers (20| , a 
negative chemical potential is necessary in order to turn 
the phantom Dark Energy a physically real hypothesis 
from the thermodynamical point of view. 

The structure of the paper is as follows. Section 2 re- 
views the basic theory of two interacting fluids with dif- 
ferent temperatures. In section 3 we apply this method 
to the case of Dark Matter and Dark Energy interaction. 
In section 4, we discuss briefly the Le Chatelier-Braun 
principle, which is related to the thermodynamical equi- 
librium of two fluids. In section 5 we analyze a simple 
interacting DE-DM model, and finally we conclude on 
section 6. 
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II. THERMODYNAMICS OF TWO 
INTERACTING FLUIDS 

In the following we review the theory of two interact- 
ing fluids, as developed by Zimdahl 2l|. The energy- 
momentum tensor of two perfect fluids (denoted by 1 
and 2) is 



dk 



with 



Ti'ik 
A 



(pa + paWu" + pAg'" , (c = 1) 



(1) 



(2) 



where pA is the energy density and pA the equilibrium 
pressure of the species A = 1,2. The 4- velocity is 
assumed to be the same for both fluids. The particle 



flow vector N\ is defined as 



UAU 



(3) 



where is the particle number density. The balance 
equation for the particle number assumes the form 



N' 



A-i 



riA + QnA = uaT / 



(4) 



where Ta is the rate of change of the number of parti- 
cles A and O = u\ is the fluid expansion rate. Particle 
production is characterized by Ta > 0, particle decay 
by < and for Ta = we have separated particle 
number conservation. 

Considering interaction between the fluids we have 
that the total energy-momentum tensor is conserved, but 
not each system separately. For a general source term 
in the energy- momentum tensor, we have 



Tiik 
A- 



A:k 



t 



A ' 



implying in the conservation equations 
PA + Q{pA + pa) = uJa 



(pa +Pa)u°' +PA.kh 



ak 



(5) 



(6) 



(7) 



where h^^ = g'^ + u^u^ . Writing the entropy per parti- 
cle as SA for each species, the Gibbs equation must be 
satisfied 



TAdsA = d— +pAd-^ 



UA 



UA 



(8) 



so that the entropy per particle variation rate is 
Uat\ {pa+Pa), 



SA 



-Ta 



tiaTa tiaTa 
The entropy flow vector S\ is defined by 



(9) 



(10) 



and we have 



^A-a 



nASA^A + nASA 
PA +PA 



SA 



nATA 



nA^A 



(11) 



where we have used relation 

The condition of energy-momentum conservation. 



0, for the system as a whole implies 



''2 J 



but there is no corresponding condition for the particle 
number as a whole, because the total number of parti- 
cles can not be conserved. We have for the total particle 
number density 



and 



where 



n = ni + n2 



h + Qn = nT . 



fiT ~ riiTi + 712 r2 : 



(12) 



(13) 



(14) 



and r is the rate by which the total particle number 
changes. 

The entropy per particle is given by 



SA 



PA +PA _ PA_ 
nATA Ta ' 



(15) 



where ^a is the chemical potential of the species A. Sub- 
stituting in Eq. pT|) yields 



The total entropy production density is 



(16) 



Qa Qa 1 Qa 

'-'■a — '^l;a + '^2;a 



M2 



Jl J2 



T 



To 



(17) 



The equilibrium condition S*.*^ = requires /ii = fi2, 
Ti = Ta and q = -q. 
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III. INTERACTION BETWEEN DARK 
MATTER AND DARK ENERGY 

Now let us apply the aforementioned results to the 
specific case of DM and DE interaction in the context of a 
FRW expansion. The Dark Energy density is represented 
by Px and is supposed to satisfy the equation of state 
(EoS) 



Px = l^Px 



(18) 



The Dark Matter energy density is represented by pdn 
and its EoS can be approximately written as 22 1 



Pdm — ndmJ^I + -ndmTdm , 



dm 



(19) 



where we have used fcs = 1. In the limit Tdm ^/ these 
reduce to pdm — ndmM and pdm — 0. Defining the inter- 
action term Uaf^rn = energy density conservation 
(|6|) of each component turns 



Pdm + '^Hpdm — Q , 

Px + 'iH{l + uj)px = -Q ■ 



(20) 
(21) 



For Q > we have the DE decaying into DM and for 
Q < is the DM that decays into DE. The cosmological 
tests indicate as many positive values |1C| as negative 
values for Q. 

Now let us analyze the process from a thermodynam- 
ical point of view. The total entropy production (|17p 
reads 



qa _ _ I Pdm p px^ -p 



Tdn 



Tx 



(22) 

Considering u) constant along the evolution, the tem- 
perature evolution for each component in terms of the 
scale factor is 



Td. 



oc a 



Tx (X a 



(23) 



thus, one expects that the present temperature of the 
Dark Energy is much greater than the temperature of the 
Dark Matter (Tx ^ Tdm)- Setting the chemical poten- 
tials as null {pdm = px ^ 0), we see that the positivity 
of the entropy production (5.'^ > 0) implies Q > 0, as 
already discussed in |12|. 

But now let us consider that the chemical potentials 
are not null. For this case the condition of positivity of 
the entropy production becomes 

1 



^ ^ / Ma: I Pd7n ^ 

Q > I TfT^x^ X + Tf; ridm^ dn 

\ X ^ dm 



{1/Tdm - l/Tx) 



(24) 



It is clear that negative values of Q are also possible, 
depending only on the signs of the chemical potentials. 
For example, if one of the components has a negative 
chemical potential while the other is null, we can have 
negative values for Q. Recently, the thermodynamical 
and statistical properties of conserved Dark Energy flu- 
ids were reexamined by considering a non-zero chemical 
potential {px 0). It was found that the entropy con- 
dition, 5*.° > 0, implies that the possible values of cu 
are heavily dependent on the value, as well as on the 
sign of the chemical potential |20[. For px > 0, the uj- 
parameter would be strictly greater than —1 (vacuum is 
forbidden) while for px < not only the vacuum but 
even a phantom-like behavior (cu < —1) is allowed. In 
any case, the ratio between the chemical potential and 
temperature remains constant for a conserved DE fluid, 
that is, Px/Tx = Pox/Tqx, where pox, Tqx are the present 
day values of the chemical potential and temperature, 
respectively. 

With basis on this discussion, taking for instance 
Mrfm = ^n.d Px = —otTx, with a a positive constant, 
the condition ([24|) turns 



Q > -aUx^xT, 



dm ■ 



(25) 



This opens the possibility of negative values of Q, so 
that the decay from DM to DE is thermodynamically 
consistent. Observe that, if DE fluid is being created, the 
quantity UxTx > 0, according to It is also interesting 
to note that this relation depends on the temperature of 
the Dark Matter, which is expected to be presently a very 
tiny value. 

As one may see, the expression p4|l is symmetric on 
the indexes x and dm . Thus the same analysis is valid 
for the case pdm 7^ and px ~ 0. Indeed, a more gen- 
eral condition necessary to have negative values to Q, 
according to is 

(26) 



Px p Pdm p 

TT^^x^ X ^ Ttt ^dm^ dn 
^ X ^ dm 



Remembering that UxTx > and Udm^dm < when DM 
is decaying into DE, for /ij, = we must have pdm > 0. 

Finally, a negative value of Q is also allowed if we re- 
quire that the entropy per particle must be conserved for 
each fluid separately. According to ([9|) we have 

Q = {Pdm + Pdm)Ydm — Mudm^dm , (27) 

Q = -{px+Px)Tx. (28) 
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and the condition turns 

It is interesting to note that, in this case, the condi- 
tion of DM decay does not depend on the decay rates, 
depending only of thcrmodynamical properties of each 
fluid. 



to the DE sector, favoring the decay in the direction of 
DE, and again implying Q < 0. Note that in this dis- 
cussion we have not imposed Q < 0. It follows from the 
action of the chemical potential, driving the direction in 
which the decay must occur. 

V. APPLICATION: Q oc Hp^ 



IV. LE CHATELIER-BRAUN PRINCIPLE 



It is also important to understand the aforementioned matter evolution law 



An interesting model of dark energy interacting 
with dark matter which has been recently proposed 
parametrizes the effect of the interaction on the dark 
13|, with Q oc Hpdm- However, 



process in the light of the Le Chatelier-Braun principle 
[2^, as recently discussed by Pavon and Wang for 
the case of DE decaying into DM {Q > 0), both with 
null chemical potentials. The principle says that when 
a system is perturbed out of its equilibrium state it re- 
acts in such manner to restore the equilibrium state or to 
achieve a new one. In the approach of null chemical po- 
tential discussed in I4I, today the answer of the system 



to the equilibrium loss is a continuous transfer of energy 
from DE to DM, corresponding to a flux / from the for- 
mer to the last one. This happens because if Q > 0, 
will increase more slowly as the Universe expands than in 
the absence of interaction and correspondingly Tdm will 
also decrease more slowly. This implies that the temper- 
ature difference between the fluids would decrease and 
eventually they will reach equilibrium. Indeed the equi- 
librium is never achieved because the expansion of the 
Universe acts as an external agent in the opposite direc- 
tion. 

Now let us see how this mechanism acts when the 
chemical potential is not null. Taking firstly the case 
with fix < and fidm — 0, we see that the negative 
chemical potential of the DE acts as a potential well, "at- 
tracting" or "sucking" the DM particles to the DE sector. 
This represents a flux F from DM to DE, in analogy to 
the previous case, but in opposite direction. Thus, if 
\F\ > I /I, the net flux is in the DE direction, represent- 
ing the decay of Dark Matter into Dark Energy, without 
violating the Le Chatelier-Braun principle. The chemi- 
cal potential acts as an "attractor" to the DM particles, 
or as an external agent driving the decay and implying 
in Q < 0. The same happens in the case fix = and 
lidm > 0. Now is the Dark Matter chemical potential 
that acts as a potential hill, "repelling" the DM particles 



although well succeeded in the background, this model 
is plagued with instabilities on the perturbations of the 
dark energy component 0|. As already shown on [2^, 
though, these instabilities are not a fairly general feature 
of these interacting models. They have shown that, in 
fact, models with Q oc Hp^ can be free of such instabili- 
ties. We rely on the conclusions of 25|, which, although 
being given for the DE decay, can be shown to be valid 
for the DM decay, too, at least for this speciflc interacting 
model. 

By assuming that baryons are conserved, the energy 
conservation law for the two interacting components 
{uaT°'l^:p = 0) is given by Eqs. ^ and At this 

point, we should mention that, although in this model 
the baryon contribution is not neglected, as usual on the 
literature, it shall not change the thcrmodynamical anal- 
ysis discussed above, as baryons are separately conserved, 
thus it has Ff, = 0, so it docs not contribute to the in- 
equality (|24|). Thus, the inclusion of baryons still allows 
for the decay of DM into DE, provided one of them has 
non-vanishing chemical potential. On another hand, the 
inclusion of baryons shall not spoil the stability of the 
perturbations, as the critical point for it is the interac- 
tion, of which baryons do not play a role. 

Thus, we shall work with the following interaction term 

Q: 



Q = 3eHpx ■ 



(30) 



where we have introduced the factor 3 for convenience. 
From this equation we can see that the sign of e deter- 
mines the sign of Q, at least while Px > 0. This model 
has already been studied on 



12 



26| 



With the interaction term ([50)) it is possible to inte- 
grate Eq. (HU, for a constant uj, and show that the 
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energy density of the DM and DE components are given 
by 

_3 , £p.oa-^[l-a-3(^+")] 

Pdm = Pdmoa H y^^) 



e + u} 



Px — PxOO- 



-3(l+e+w) 



(32) 



where we can see that in the hmit of vanishing e, we 
recover the standard, non-interacting, behaviors. As 
baryons are conserved, we have pt = pho(l + z)'^, so the 
total pressureless matter density pm (= Pb + Pdm) goes 
like: 



Pm — PmOO, 



(33) 



We choose to work on a spatiaUy flat Universe, which 
is predicted by inflation and is in agreement with the 
(Wilkinson Microwave Anisotropy Probe) WMAP obser- 
vations 



27 



28 



29| . Thus, neglecting the radiation contri- 



bution, the Friedmann equation for this modified XCDM 
cosmology can be written as: 



H 



LO - 



■(1 + ^) 



(l + z)3 + 



3(l+e+i^) 



(34) 



where i?, fif, and Vim (= l^ — ^x) are, respectively, the 
Hubble parameter, the baryons and matter present den- 
sity parameters and we have dropped the subscript '0' 
for convenience. 

This model differs from the model of 



12 



26[, as they 



neglect the baryon contribution. Our aim here is to use 
the most recent observations and the inclusion of baryons 
in order to put stringent constraints on this model. 

We fix fif, = 0.0437 from the WMAP 5-year observa- 
tions |29| , which is also in agreement with the nucleosyn- 



thesis predictions, so we have four free parameters on our 
model, namely, h, Q.m (= f^dm + f^b), ^ and e (neglecting 
the radiation contribution, which is of order Vl^h? ^ 10^^ 



28|). The dark energy density parameter Vlx is not a 
free parameter as we choose to work on a spatially flat 
Universe. We shall combine various cosmological obser- 
vations in order to constrain these parameters. 

The Baryon Acoustic Oscillations (BAO) given by the 
acoustic oscillations of baryons in the primordial plasma, 
leave a signature on the correlation function of galaxies as 
observed by Eisenstein et al. (2005) 30[. This signature 
furnishes a standard rule which can be used to constrain 



the following quantity [30 1 



A 



1/2 



n1/3 



-r(z,) 



2/3 



0.469 ±0.017, (35) 



where H = H[z)/Hq, which is given by Eq. ([34|) . z, = 
0.35 is a typical redshift of the SDSS sample, and r(z*) 
is the dimensionless comoving distance to the redshift z^, . 
This quantity has been constrained by Eisenstein et al. 



30j to be ^ = 0.469 ± 0.017 and we shall use a gaussian 



prior on it. 

SNe la luminosity distances constrain mainly the dark 
energy density and dark energy EoS, and we choose to use 
the up-to-date greatest combination of known SNs data 
sets, namely the Union compilation 3l|. This sample 
is composed of 307 SNs, which includes the recent large 
samples of SNe la from the Supernova Legacy Survey 
3 and ESSENCE Survey 3^], older datasets, as weU 
as the recently extended dataset of distant supernovae 
observed with Hubble Space Telescope (HST). We choose 
to use the "without systematics" Union dataset. We use 
a fairly well known [33[ analytical marginalization over h 
to combine the SNe luminosity distances data. 

An useful quantity to characterize the position of the 
Cosmic Microwave Background (CMB) power spectrum 
flrst peak is the shift parameter, which is given, for a flat 
Universe, by [34| : 

(36) 



TZ = \ nl 

^(^) 

where Zr is the recombination redshift. The recombina- 
tion redshift can be calculated by using the fitting for- 
mulas by Hu and Sugiyama [35[. The CMB shift pa- 
rameter given by the WMAP 5-yr Monte Carlo Markov 
Chain (MCMC) analysis, assuming a standard (non- 
interacting) XCDM model is 7e = 1.710 ± 0.019. 

As we marginalize over h on the SNe data analysis, 
the only dependence on the Hubble constant comes from 
the shift parameter, via recombination redshift. But this 
dependence is too weak, so we choose to fix h from the 
HST Key Project best fit, h = 0.72±0.08 0. So, finally, 
we have to minimize the following quantity: 

2 , _-, „s 2 



X 



A - 0.469 
0.017 



7^-- 1.710 
0.019 



(37) 



where Xsn ^^"^ marginalized from the SN data, 
given by Xgp^ ~ —2 In Cswdh, Csn is the SNe mag- 



nitude likelihood, Csn = exp 



_ 1 Y^3U7 
2 2^i=l 



6 




FIG. 1: a) Plane Q^n. - £ showing the 68.3%, 95.4% and 99.7% contours of probabihty, as well as the £ = line (dashed line), 
b) Marginalized likelihood, normalized at the peak (solid line) and CDF (long-dashed line) for the e parameter, as well as the 
e = line (short-dashed line) and P{e > 0) line (dotted line). See text for details. 



fi is the distance modulus, given as a function of the lu- 
minosity distance cIl in Mpc as = 51og((ii) -I- 25. 

The total likelihood is given by £ oc e"^ Marginal- 
izing C over Lu, we can see how much e and f2„i can be 
constrained by the data. The 1, 2 and 3 — tr contour 
levels are shown on Fig. [1^. 

We have found, then, from this analysis, Qm = 



n 90/1 +0.019+0.039-1-0.059 „^ i ^ 
U.Z01_Q g2o_o. 038-0. 056 '^'^^^^ ^ 



n r)9R+0. 021+0. 038+0. 051 



for 1, 2 and 3(7, respectively. We have obtained a relative 
X^/i^ = 311.66/305 = 1.02, where v is the number of de- 
grees of freedom. As one can already see on Fig. [1^, the 
possibility of e < is clearly favored over the possibility 
of e > 0. As one may also see, the standard value, e = 0, 
is marginally inside the 1-a region, in the joint analy- 
sis. In order to give a more quantitative result, we also 
marginalize over flm with the aim to find the likelihood 
of £, C{e). This is shown on Fig. [TJd. 

As one may see, the largest area under this curve is 
given for e < 0. It can be confirmed by calculating the 
probability P{e < 0), which is given by 



Pie < e') = 



C{e)de 



(38) 



This is nothing more than the cumulative distribution 
function (CDF) of e and it is also shown on Fig. [TJj. By 
putting £r' = on (IMl), we find that P{e < 0) = 93.27%. 
It corresponds to a ratio P{e < 0)/P{e > 0) = 13.86, 
and an almost 2-cr c.l. in favor of e < in comparison 
with e > 0. Here, on the marginalized analysis, we can 



also see that the standard value, e = 0, is outside the 1-a 
region, but is well inside the 2-a region, though. 

VI. CONCLUSIONS 

Assuming that Dark Energy is amenable to a fluid de- 
scription with a well defined temperature and chemical 
potential, we have analyzed the possibility of having DM 
decaying into DE through thermodynamical considera- 
tions. We have shown that, once at least one of the fluids 
is provided with a non null chemical potential, this pos- 
sibility is allowed by the second law of thermodynamics. 
We have also showed that the Le Chatelier-Braun prin- 
ciple remains valid. Next we have applied this result to 
put observational constraints on a fairly simple model of 
interaction between DM and DE. We have found that, in 
the context of this model, the decay of DM into DE is 
reasonably favored over the inverse process. This model 
has the correct behavior of being DM dominated in the 
past and DE dominated in the future. Our result gives 
a theoretical basis and opens the possibility of extending 
and analyzing other models of interacting DE present 
on the literature. However, it should be mentioned that 
our results rely on the assumption of a nearly standard 
evolution of per tubations on interacting DE models, as 
indicated by 25[. While this result applies to the model 
pop , a more general perturbation analysis is still lacking. 
Work in this direction is being done and can appear in a 
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forthcomming communication. 
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